Fundamental and vortex solitons in a two-dimensional optical lattice 
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Fundamental and vortex solitons in a two-dimensional optically induced waveguide array are 
reported. In the strong localization regime, the fundamental soliton is largely confined to one 
lattice site, while the vortex state comprises of four fundamental modes superimposed in a square 
configuration with a phase structure that is topologically equivalent to the conventional vortex. 
However, in the weak localization regime, both the fundamental and vortex solitons spread over 
many lattice sites. We further show that fundamental and vortex solitons are stable against small 
perturbations in the strong localization regime. OCIS codes: 190.0190. 



Discrete optical spatial solitons are localized modes 
of nonlinear waveguide arrays that form when discrete 
diffraction is balanced by nonlinearity [1]. When a low 
intensity beam is focused into a waveguide array, the 
propagating field spreads over many sites (due to optical 
tunneling) exhibiting a typical discrete diffraction pat- 
tern with the intensity mainly concentrated in the outer 
lobes [2]. However, at sufficiently high power, the beam 
self-traps to form a localized state. 

Discrete solitons in waveguide arrays were first pre- 
dicted to exists as solutions to the discrete nonlinear 
Schrodinger equation [3] and later observed in AlGaAs 
waveguide arrays [2]. This experimental observation 
stimulated much new research such as diffraction man- 
agement [4,5], soliton interaction and beam steering [6], 
and discrete solitons in a two-dimensional (2D) photore- 
fractive optical lattice [7-9]. The experiments in [7-9] 
are particularly interesting since the 2D waveguide ar- 
ray there is formed optically, thus it is very versatile and 
easily tunable. This allowed a host of 2D localization 
phenomena to be observed in [7-9]. However, many ques- 
tions on discrete solitons in 2D optical lattices still remain 
open. For instance, stability properties of such solitons 
have not been carefully studied. In addition, new local- 
ization states such as those analogous to vortex solitons 
in a homogeneous waveguide have not been explored. 

In this Letter, fundamental and vortex solitons in a 
2D optical lattice are found and their stability analysed. 
In the strong localization regime, the fundamental soli- 
ton is confined largely on one lattice site with a uniform 
phase. The vortex soliton comprises of four fundamen- 
tal modes located at the bottoms of the optical potential 
in a square configuration with a phase structure that is 
topologically equivalent to the conventional vortex. By 
winding around the zero intensity position along any sim- 
ple closed curve, the phase of the vortex state acquires 
2-7T increment. We call this structure a vortex cell. When 
the localization is weak, both fundamental solitons and 
vortex-cells spread over many lattice sites. In the strong 
localization regime, we show that the fundamental soliton 
as well as the vortex cell are both linearly and nonlinearly 
stable under weak perturbations. 

We start with the 2D nonlinear Schrodinger equation 
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where V = Vq (cos 2 X + cos 2 Y) is the optical lattice po- 
tential, and Vo is its intensity. Such a potential can be 
obtained by optically interfering two pairs of laser beams 
[9]. Transverse distances X and Y in (I) have been nor- 
malized by the lattice period D, and distance z normal- 
ized by AuqD 2 /ttA, where hq is the refractive index, and 
A is the wavelength of the beam [9]. For typical values 
D = 9/iTO, no = 2.3 and A = 0.5/im [8,9], the unit dis- 
tance z in (I) corresponds to 0.5mm in physical distance. 
Without the lattice potential, solitons would suffer col- 
lapse under small perturbations [12]. However, as we 
shall show here, optical lattices can suppress the collapse 
of fundamental solitons and vortex-cells. 

Eq. (1) conserves two quantities: the power P = 
IZoI- \i>\ 2 dXdY, and the energy E: 
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We look for stationary solutions of the form ip(X, Y, z) = 
e~ lfJ,z u(X,Y) where /i being the propagation constant of 
the soliton. Then u(X,Y) satisfies 
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Solutions to this equation can be obtained by a Fourier- 
iteration method as was suggested by Petviashvili [14]. 

Fundamental solitons. A fundamental soliton of Eq. 
(3) has a single main hump sitting at a bottom of the 
potential, say (X,Y) = (ir /2,ir/2). Two examples corre- 
sponding to propagation constants /j, = and 0.88 with 
Vo = 1 are displayed in Fig. l(c,d). We see that for 
small fi (Fig. lc), the beam is largely confined on one 
lattice site, while at higher n (Fig. Id), it spreads over 
many lattice sites. To quantify these solitons, we cal- 
culate the dependence of the normalized power P on 
the propagation constant /i for Vo = 1 as displayed in 
Fig. la. When fi —> — oo, P approaches a constant 11.70. 
This is apparently because in this limit, the fundamental 
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soliton is highly localized, thus it approaches the lattice- 
free fundamental-soliton state, which has critical power 
P c w 11.70. As \i goes to a cut-off value which is approx- 
imately 0.95, P appears to go to infinity. In this limit, 
the fundamental state becomes uniformly distributed in 
space. Thus, this cut-off value should be the boundary 
of the band gap in the linear-wave spectrum [10]. When 
M « 0.72, P is minimal. When fj, > 0.72, dP/dfi > 0. In 
this region, the Vakhitov-Kolokolov (VK) theorem sug- 
gests that the soliton is linearly unstable due to the pres- 
ence of a purely-real unstable eigenvalue a in the lin- 
earized equation [11]. We have confirmed this instability 
by numerically simulating the linearized version of (1) 
around the above soliton. The results for unstable eigen- 
values are shown in Fig. lb. On the other hand, when 
fi < 0.72 (Vb = 1), the VK instability is absent, and the 
soliton state is linearly stable. 

For the two-dimensional self-focusing case, collapse is 
an important issue. The above linear stability analy- 
sis does not guarantee that the fundamental soliton will 
not collapse under small perturbations. In the study of 
collapse, the energy E plays an important role. In the 
absence of the lattice potential [12] or when the poten- 
tial is harmonic [13], the soliton collapses if its energy 
is negative. In case the energy is positive, however, the 
soliton collapses only if it is strongly perturbed. For the 
optical lattice, we have calculated the energy E of fun- 
damental solitons at various values of fj, and plotted the 
results in Fig. 1(a). The energy is found to be always 
positive. Thus we can expect that this state is able to 
withstand small perturbations without collapse. To con- 
firm the above expectations on stability, we numerically 
study the nonlinear evolution of the fundamental soliton 
under small perturbations by directly simulating Eq. (1) 
with initial condition 
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where e « 1, and u p (X,Y) is the initial perturbation. 
We first take u p to be white noise. A large number of sim- 
ulations with small e and various realizations of random- 
noise perturbations have been performed, and we have 
found that for Vb = 1, if H < 0.72, the fundamental soli- 
ton is indeed stable against white noise perturbations; 
when li > 0.72, the soliton is unstable. To study the 
nonlinear evolution process, we now take u p = 1. For 
Vo = 1, fj, = and e = ±0.01, nonlinear evolutions are 
plotted in Fig. 2a. We see that the perturbed soliton only 
oscillates weakly around the fundamental-soliton state, 
meaning that the soliton is both linearly and nonlinarly 
stable. On the other hand, at Vb = 1 and /i = 0.88 (where 
the soliton is linearly unstable), the dynamics is different 
as two scenarios are identified: (i) at higher input power 
(e > 0), the perturbed state relaxes into a z-periodic 
bound state (Fig. 2b, c); (ii) at lower input power (e < 0), 
the perturbed state decays into linear Bloch waves (Fig. 



2b, d) (similar scenarios can be found in [15] for a different 
system). Collapse is not observed in either scenarios. 
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FIG. 1. (a) Normalized power P and energy E of funda- 
mental solitons versus fx for Vo = 1. (b) unstable eigenvalues 
a of these solitons for Vo = 1; (c, d) profiles of fundamental 
solitons at fx = and 0.88 (Vb = 1) respectively. 



(a) 



a. 2 

I 1.9 





e=0.01 




e=-0.01 



(b) 

£=0.01 



e~ 0.01 



20 

t 




50 100 
t 

(d) 



-4jtl 
-2jll 



2itl 
4x1 



Hi 



-4jl-2lt 271 4ti 671 
X 



-47C-27C o 27c 4ti 6tc 
X 



FIG. 2. Nonlinear evolutions of fundamental solitons under 
perturbations (4) for Vo = 1, e = ±0.01, u p = 1 and fx — 
(a), 0.88 (b). Snapshots of the soliton intensity corresponding 
to the dynamics depicted in (b) for e = 0.01, z — 80 (c) and 
e = -0.01, z = 60 (d). 

Vortex cell. In addition to the fundamental solitons, 
we have numerically found vortex solitons as well. Two 
examples with Vb=l, fx = and 0.82 are shown in Fig. 3 
(b,c,d). At fx = (strong localization regime), the vor- 
tex state comprises of four fundamental solitons super- 
imposed in a square configuration with a phase structure 
which is topologically equivalent to a conventional vortex 
[see Fig. 3b, c]. By winding around the center along any 
closed curve, the phase of the vortex acquires 27r incre- 
ment, thus, we name it a vortex cell. At fx = 0.82 (weak 
localization regime), the vortex cell spreads out to more 
lattice sites and becomes more intricate, as can be seen 
in Fig. 3d. But its phase structure is almost the same 
as with fx = 0. We should point out that these vortex 
cells are different from conventional vortices without op- 
tical lattice on a major aspect: their densities and phases 
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depend on both r and 9. The normalized power and en- 
ergy diagrams versus propagation constant for these vor- 
tex cells at Vo = 1 are shown in Fig. 3(a). This figure 
is similar to that for fundamental solitons [Fig. 1(a)], 
except that both P and E here are about four times 
larger. Unstable eigenvalues a of vortex-cells are deter- 
mined by simulating the linearized equation (1) around 
vortex-cells. The results are shown in Fig. 4a for Vo = 1. 
We see that vortex cells experience an oscillatory instabil- 
ity for /i > —7.8 and become stable for /i < —7.8. These 
vortex cells also suffer the VK instability in the region 
fj, > 0.73 where dP/d^i > [see Fig. 3a]. However, the 
oscillatory instability is much stronger as it occurs over 
a wider region and has a higher growth rate. Fig. 3(a) 
shows that vortex cells also have positive energy. Thus, if 
a vortex cell is linearly stable, it should be able to resist 
collapse under small perturbations [13]. However, if it 
suffers the linear oscillatory instability discussed above, 
this instability could result in power exchange from one 
part of the cell to another so that the intensity at some 
small spots becomes high, triggering local collapse to oc- 
cur. We have observed this scenario numerically. An 
example is shown in Fig. 4(b,c,d) which displays the de- 
velopment of a vortex cell with /j, = and Vb = 1 when 
it is amplified by 1% initially. 



stabilize both fundamental and vortex solitons. 



(a) (b) 




-2k -% o TC 2k -2k -71 K 2k 
X X 



FIG. 4. (a): Unstable eigenvalues of vortex cells with 
Vo — 1; (b,c,d) instability development of the vortex with 
Vo — 1 and n — [see Fig. 3(b,c)] when it is initially ampli- 
fied by 1%; intensity plots at z = 50, 51.5 and 52 are shown 
respectively. 

In conclusion, we have studied new types of funda- 
mental and vortex solitons in a two-dimensional optical 
lattice potential, and shown that both solitons are stable 
in the strong localization regime. 

The work of J.Y. was supported by NSF and NASA. 
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FIG. 3. (a) Normalized power P and energy E of vortex 
cells versus /i for Vb = 1. (b,d) intensity plots of vortex cells 
with Vo = 1, fj, = and 0.82 respectively, (c) the phase plot 
of the vortex cell in (b). 

We next discuss the effect of varying the potential 
strength Vo on the formation and stability of fundamen- 
tal and vortex solitons. For this purpose, we have chosen 
Uo = 1.5, and repeated most of the above calculations. 
The results are summarized as follows: (i) fundamental 
solitons and vortex-cells exist at higher values of /j, (up 
to approximately 1.35); (ii) at the same value of /i, both 
the fundamental and vortex solitons have lower power 
(P) and higher values of energy (E); (hi) the oscillatory 
instability suffered by vortex cells is reduced (it could 
even be completely suppressed when U becomes even 
larger). Thus, we conclude that higher lattice potentials 
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